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1 Introduction 



Recently it has been shown — as it had been already realized in old times [1] — that not 
only nonlinear, but also a large class of linear equations (including, in particular, the 
wave equations) admit of "soliton-like" solutions. Those solutions [2] are localized, and 
travel along their propagation axis practically without diffracting (at least until a certain 
field-depth [2, 3, 4]): Such wavelets were indeed called "undistorted progressing waves" by 
Courant and Hilbertfl]. Let us recall that their peak- velocity V can assume any val- 
ues[5,6,2] < V < oo, even if we are mainly interested here in their localization properties 
rather than in their group-velocity. In the case of wave equations, the localized solutions 
more easy to be constructed in exact form resulted to be the so-called "(superluminal) 
X-shaped" ones (see Refs. [4,7,8,2], and refs. therein). 

The X-shaped waves, long ago predicted [6] to exist within Special Relativity (SR), 
have been first mathematically constructed[9,2] as solutions to the wave equations in 
Acoustics [4], and later on in Electromagnetism (namely, to the Maxwell equations[7]), 
and soon after produced experimentally [10]. Only very recently, subluminal localized 
solutions have been suitably worked out in exact formfll], even for the case of zero speed 
("Frozen Waves"). [12] 

It was soon thought that, since the mentioned solutions to the wave equations are non- 
diffractive and particle-like, they may well be related to elementary particles (and to their 
wave nature) [13,14]. And, in fact, localized solutions have been found for Klein-Gordon 
and for Dirac equations[13,14]. 

However, little work[15] has been done, as far as we know, for the {different) case 
of the Schroedinger equation^ Indeed, the relation between the energy E and the im- 
pulse magnitude p = |p| is quadratic [E = p 2 /(2m)] in the non-relativistic case, like 
in Schroedinger's, at variance with the relativistic one. But, as we were saying, the 
nondiffracting solutions, which are essentially superpositions of Bessel beams and are 
currently called Localized Waves, would be quite apt at describing elementary particles: 
much more than the gaussian waves. In this paper we show that indeed, mutatis mutan- 
dis, Localized Solutions exist even for the ordinary Schroedinger equation within standard 
Quantum Mechanics; and we obtain both approximate and exact solutions, also setting 

*For some work in connection with the ordinary Schroedinger equation, see for instance, besides [7], 
also Refs. [14]. 
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forth for them particular examples. In the ideal case such solutions bear infinite energy, 
as well as spherical or plane waves: we shall therefore show how to obtain finite-energy 
solutions. At last, we shall briefly consider solutions for a particle moving in the presence 
of a potential. 

Before going on, let us recall that, in the time-independent realm — or, rather, when 
the dependence on time is only harmonic, i.e., for monochromatic solutions — , the (quan- 
tum, non-relativistic) Schroedinger equation is mathematically identical to the (classical, 
relativistic) Helmholtz equation[16]. And many trains of localized X-shaped pulses have 
been found, as superpositions of solutions to the Helmholtz equation, which propagate, for 
instance, along cylindrical or co-axial waveguides [17]; but we shall skip all the cases[18] 
of this type, even if interesting, since we are concerned here with propagation in free 
space, even when in the presence of an ordinary potential. Let us also mention that, in 
the general time- dependent case, that is, in the case of pulses, the Schroedinger and the 
ordinary wave equation are no longer mathematically identical, since the time derivative 
results to be of the fist order in the former and of the second order in the latter. [It has 
been shown that, nevertheless, at least in some cases [19], they still share various classes 
of analogous solutions, differing only in their spreading properties [19]]. Moreover, the 
Schroedinger equation implies the existence of an intrinsic dispersion relation even for 
free particles. 

Another difference, to be kept here in mind, between the wave and the Schroedinger 
equations is that the solutions to the wave equation suffer only diffraction (and no dis- 
persion) in the vacuum, while those of the Schroedinger equation suffer also (an intrinsic) 
dispersion even in the vacuum. 

Let us repeat that the majority of the ideal localized solutions we are going to con- 
struct are endowed with infinite energy. We shall treat also a finite-energy casa[j] only 
towards the end of this paper: In fact, infinite-energy solutions themselves, even without 
truncating them in space and time, results to be rather useful for describing wavepackets 
in regions not too extended in the transverse direction; as we shall see below. 



'In such cases the solutions travel undistorted and with a constant speed along a finite depth of field 
only. 
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2 Bessel beams as localized solutions (LS) to the 
Schroedinger equation 

Let us consider the Schroedinger equation for a free particle (an electron, for example) 

w + ?££ = o. (., 

If we confine ourselves to solutions of the type 

their spatial part F obeys the reduced equation 

V 2 F + k 2 F = , (2) 

with k 2 = p 2 /h 2 and p 2 = 2mE (quantity p = |p| being the particle momentum, and 
therefore k = |fc| the total wavenumber). Equation (2) is nothing but the Helmholtz 
equation, for which various simple localized-beam solutions (LS) are already known: In 
particular, the so-called Bessel beams [2], which have been experimentally produced since 
long [20]. 

Namely, let us now look — as usual — for solutions (cylindrically symmetric with 
respect to [w.r.t.] the z-axis) of the form 

ijj(p,z;t) = R(p) Z(z) T(t) , 

and explicitly indicate, mainly for clarity's sake, the subsequent steps. Equation Q then 
becomes 



J^c^ ( dR\ l d 2 Z 
pR dp \ dp ) z dz 2 



2mldT 



so that 

i2m dT 



HT dt 



-k 2 T = e- lEt/h (4) 
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where, let us repeat, E = p 2 /(2m) = k 2 h 2 /(2m) (and in fact the last exponential is often 
written as exp [— icut]). 



Analogously, we have 



J^8_ ( 8R\ 2 ld 2 z 

pRdp\ P dp) +k ~ zdz 2 {b) 



and therefore 



-ig = k 2 z =► * = e**-, (6) 

where the constant k z = fcy = = is the longitudinal wavenumberj*] We will 
suppose k z > 0, that is, p z > 0, to ensure that we deal with forward traveling beams only. 
As a consequence, the (transverse) function R = R(p) obeys the equation 

\h ("f ) + (p " = (7) 

which is a Bessel differential equation admitting as solution the Bessel function^] 

R(p) = J (pk p ) , (8) 
where the constant k p = k± = p±/h is the transverse wavenumber, and 

k 2 p = k 2 -k 2 = 2mE/h 2 - k 2 z . (9) 

To avoid any divergencies, it must be k 2 > 0, that is, k 2 > k 2 ; namely, it must hold [see 
(a) in Figjl] the constraint 

2 

E > — (10) 
~ 2m K J 

[Notice that, to avoid the appearance of evanescent waves, one should postulate k z to be 
real; but such a condition is already included in our previous assumption that k z > 0]. 
In the following, to simplify the notations, we shall also put [p = hk ] : 



■''Since the present formalism is used both in quantum mechanics and in electromagnet ism, with a 
difference in the customary nomenclature, for clarity's sake let us here stress, or repeat, that k = p/Ti; 
k p = k± = p±/h; cu = E/h; while k z = fcii = P\\/Ti = p z /h is often represented by the (for us) ambiguous 
symbol j3. 

§The other Bessel functions are not acceptable here, because of their divergence at p = or for p — > oo. 
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The solution is therefore: 



P P = Pi 



ip{p,z;t) = J (pp p /h) exp[i(zp z - Et)/h] 



(11) 



together with condition (9). Equation (11) can be regarded as a Bessel beam solution 
to the Schroedinger equation. This result is not surprising, since — once we suppose the 
whole time variation to be expressed by the function exp [iut] — both the ordinary wave 
equation and the Schroedinger equation transform into the Helmholtz equation. Actually, 
the only difference between the Bessel beam solutions to the wave and to the Schroedinger 
equation consists in the different relationships among frequency, longitudinal, and trans- 
verse wavenumber; in other words (with E = uh): 



In the case of beams, the experimental production of LSs to the Schroedinger equa- 
tion can be similar to the one exploited for the LSs to the wave equations (e.g., in Optics, 
or Acoustics): Cf., e.g. Figure 1.2 in the first one of Refs.[8], and refs. therein, where 
the simple case of a source consisting in an array of circular slits, or rings, were consid- 
ered]^] In the Table we refer to a Bessel beam of photons, and a Bessel beam of (e.g.) 
electrons, respectively. We list therein the relevant quantities having a role, e.g., in Elec- 
tromagnetism, and the corresponding ones for the Schroedinger equation's spatial part 
h 2 'V 2 F + 2m E F = 0, with F = R(p) Z(z) . The second and the fourth lines have been 
written down for the simple Durnin et al.'s case, when the Bessel beam is produced by 
an annular slit (illuminated by a plane wave) located in the focus of a lens [20]. 



^For pulses, however, the generation technique must deviate from Optics', since in the Schroedinger 
equation case the phase of the Bessel beams produced through an annular slit would depend on the 
energy. 



p 2 p = E 2 /c 2 —p 2 z for the wave equation; 
p 2 = 2mE — p 2 for the Schroedinger equation. 



(12b) 
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WAVE EQUATION 


SCHROEDINGER EQUATION 


k — - 

C 

IC <~**~ ' — X? 

L.2 . .. ui 2 _ p 


p = sj2mE 
Pp-jP 
p 2 p = 2mE — p\ 
p 2 z = 2mE{\ - £) 



In this Table, quantity / is the focal distance of the lens (for instance, an ordinary lens 
in optics; and a magnetic lens in the case of Schroedinger charged wavepackets), and r 
is the radius of the considered ring. [In connection with the last line of the Table, let 
us recall that in the wave equation case the phase- velocity u/k z is almost independent 
of the frequency (at least for limited frequency intervals, like in optics), and one gets a 
constant group- velocity and an easy way to build up X-shaped waves. By contrast, in the 
Schroedinger case, the phase-velocity of each (monochromatic) Bessel beam depends on 
the frequency, and this makes it difficult to generate an "X-wave" (i.e., a wave depending 
on z and t only via the quantity z — Vt) by using simple methods, as Durnin et al.'s, based 
on Bessel beams superposition. In the case of charged particles, one should compensate 
such a velocity variation by suitably modifying the focal distance / of the Durnin's lens, 
e.g. on having recourse to an additional magnetic, or electric, lens.] 

Before going on, let us stress that one could easily eliminate the restriction of axial 
symmetry: In such a case, in fact, solution (11) would become 

1>(p,z,<p;t) = J n (pp P /h)e izp °/ n e- iEt / n e in * , 

with n an integer. The investigation of not cylindrically-symmetric solutions is interesting 
especially in the case of localized pulses (cf. Sect. 3): and we shall deal with them below. 

3 Localized pulses as solutions to the Schroedinger 
equation (approximate method) 

Localized (non-dispersive, besides non-diffracting) pulses can be constructed, as solutions 
to the Schroedinger equation, both by having recourse to the standard "paraxial approx- 
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Figure 1: The parabola and the chosen straight-line have equations E = pi /(2m) and 
E = Vp z , respectively. The intersection of our straight-line with the parabola corresponds 
to the value E = 2mV 2 . The allowed region is the one internal to the parabola, since it 
must be E > pi/ (2m). 

imation", and in an exact, analytic way. Let us start with the approximate method. 

Let us go back, then, to our Bessel beam solution (11), with condition (10). We 
can obtain localized (non-dispersive) pulses, as solutions to Schroedinger's equation, by 
suitably superposing the beam solutions (11), and by selecting in the plane (p z ,E) the 
straight-line [see Fig.l]: 



vith V a chosen constant speed; so that from eq.(10) one gets the important condition 



E 



V Pz ; 



(Pz > 0) 



(13) 



E < 2mV 2 



(14) 



and eq.(ll) can consequently be written 



(ir) 



where now p 2 p = (2mE — p 2 z ) = E(2m — E/V 2 ) and we introduced the new variable 



c 



vt . 



(15) 
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Localized-wave solutions can be therefore obtained through the superposition (see 



the weight-function S(E) being a suitable energy-spectrum (with the dimensions, as usual, 
of the inverse of an Energy), while Af is a "normalization" constant which normalizes 
to 1 the peak- value of 2 and (since it multiplies a dimensionless integral) bears the 



noted that we are integrating, in the space (p z ,E) along the straight-line (13), that is, 
E = Vp z . This corresponds to superposing Bessel beams all endowed with the same 
phase- velocity Vph = V. The resulting pulse will possess V as its group-velocity (namely, 
as its peak- velocity) , since it is well-known that, when the phase- velocity V^h does not 
depend on the energy or frequency, the resulting pulse happens to travel with the group- 
velocity = duj/dk z = V p h = V: cf. refs. [17,2,21] and refs. therein. Due to constraint 
(14), we are actually integrating along our straight-line from to 2mV 2 (see Fig.l). 

It is important also to note explicitly that each solution C) given by eq.(16), 
depending on z (and t) only via the variable ( = z — Vt, does represent a pulse that appear 
with a constant shape to an observer traveling with speed V along the wave motion-line z: 
in other words, it represents a pulse which propagates rigidly along z. Therefore, eqs. (16) 
are already — as desired- non- dispersing and non- diffracting ("localized") solutions to the 
Schroedinger equation. 

Integrals (16), however, appear difficult to be analytically performed, independently 
of the spectrum S(E) chosen. 

To overcome this difficulty, let us rewrite eq.(ll') as a function of p p only, by exploiting 
eq.(12b), which can be written E 2 /V 2 — 2mE + p 2 p = 0, and yields 



Fig.l): 




(16) 



dimensions [J\f] = [L 3 ^ 2 ], to respect the ordinary meaning 

of |tf(p,C)| • It should be 




(17) 



where 



Ppmax = mV , 
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as it comes by deriving eq.(12b) with respect to E. 



Therefore, eq.(ll') becomes 

*P(P,0 = Jo( PPp /h) exp[z^Cv / i r J^ ] S( Pp /n) e^C (11") 

with < p p < p pmax , wherejjijlet us repeat, p pmax = mV . Then, the Localized Solutions 
will be written as 



¥(p,C) = Me imV ^ h J Q mV dp p U P -^)S{p p ) exp 



imV I _ p 2 p 



h V m?V 2 



Let us notice that, in the new variable p p , the Bessel function, previously written as in 
eq.(16), gets, as we have seen, the simplified expression Jo(pp P )- 

It is now enough to choose a weight-function S that is strongly bumped around the 
value p p , in the interval [0 , m V^] , with 

p p < mV , (19) 

for being able to integrate from to oo with a negligible error. Namely, let us now adopt 
the so-called paraxial approximation. Under condition (19), one can approximate the 
exponential factor as follows: 

p 2 1 p 2 

mV\ 1 -t— ~ mV — 



m 2 V 2 2 mV ' 

so that eq.(18) can be eventually written in terms of an integration from to oo: 

2 

tf(p, C) = Afe 2imV « n | o °° dp p Jo(^) S( Pp ) exp [i^yC] ■ (20) 

Let us now examine various special cases of weight-functions S(p p ) obeying the pre- 
vious conditions: that is, well localized around a value p p mV. 



"For the sake of clarity, let us repeat that, when the phase- velocity V becomes (as in our case) the 
group-velocity, V g = V , then the component p p of p acquires mV as its maximum value. It holds, 



moreover, Jp 2 — p 2 p = p\\ = p z , which just equals p, since in the present case V = \V\ = V z . 
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3.1 Some examples of approximate Localized Solutions to the 
Schroedinger equation (paraxial approximation) 

As already claimed, we are for the moment adopting the paraxial approximation, since it 
yields good, and interesting enough, results: Only in the subsequent Sections we shall go 
on to the exact, analytical approach. 

First of all, let us consider the simple spectrum 

S(P P ) = ±qp P e- qp2 » (21) 
(with the dimensions, now, of the inverse of an Impulse), with 

q = -^77- 2 (22a) 

so that the above conditions merely imply the dimensionless constant a to be 

a > 1 . (22b) 

In this case, also the total spectral-width Ap p results to be Ap p <C mV: and this too 
supports the fact that our integral can indeed run till oo. In eq.(20), one can then perform 
(analytically) the integration, and get the solutions 

*(p,C) * M 4# e—</» ± exp h 4%fi ^ ] , (23) 
quantity q being still the one defined in eq.(22a), with a >> 1; while function Q is 

Equation (23) constitutes an interesting solution of the Schroedinger equation: It describes 
a wavepacket rigidly moving with the chosen speed V. The maximum of its intensity 
|\l/| 2 occurs at 

P = 0; C = , 
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and therefore also such a maximum travels with the speed V, as expected (since ( = 
z - Vt). For C = one gets [a > 1]: 

|*(p,C = 0)| 2 ^AfU exph^j], (25) 
and the transverse localization Ap of the wavepacket results to be 

Ap=A^, (25') 

which shows also the role of a (and therefore of q) in regulating the wavepacket (constant) 
transverse total width. 

By contrast, putting p = into eq.(23), we end up with the expression [still with 
a > 1]: 

I*(P = 0,C)| 2 ^4 2 f^_ , (26) 

which corresponds to 

a > r~n — 7 2a/i 
AC = Ve 2 - 1 



mV 

Solution (26) is represented in Fig. 2. 



Let us briefly consider a few further possible spectra. We shall go on confining 
ourselves to the simple case of cylindrical symmetry, but analogous solutions can be 
easily found also for more general non-symmetrical cases. 

As the second option, let us choose the new spectrum 

S(p p ) = - e-«* , (27) 
Pp 

quantity q being defined in eq.(22a), and condition (22b) being enforced, so that q ^> 
l/(m 2 V 2 ) and, again, Ap p <^ mV . Equation (20) yields the new solution 
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Figure 2: Behavior of = 0, C) I m ec L-(26), as a function of (/(2hqmV). 

C)^^-7(0, ^ J exp [— C] , (28) 

where function Q is defined in eq.(24), and 7, here, is the "incomplete gamma func- 
tion".^] 

7(0,4) = - 1 (-1,A)-A~ 1 e- A 

with 

j(-l,A) = -A- l e~ A $(1,0; 4) 

= -A- l e- A [l-§(l^A)} , 

function $ being the "Probability Integral" , that in the present case can be defined as 

1 f°° , a — e~ Ax 



4(1,0; -4) = m r°dx 



r(i) y 1-e-* 

The maximum, also for solution (27), occurs at p = ( = 0. 
v4s a third option, we choose 
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S(P P ) = q Pp e-^I ( S -f: 



(29) 



always with a 1, quantity q being given by eq.(22a), s a constant with the dimensions of 
a Length (regulating the spectrum bandwidth), and J being the Modified Bessel Function; 
one gets from eq.(20) the further new solution 



qh a 
2Q 



s 2 -p 2 
4Q 



Jo 



sp_ 
2Q 



(30) 



As the last option, let us choose 



S(p P ) = qp P e qPp M S P P ) , 
from eq.(20) it follows the fourth solution 



(31) 



*(p,C) ^ M ^e » c exp 



s 2 + p 2 
4Q 



20, 



(32) 



4 Exact Localized Solutions to the Schroedinger 
equation (for arbitrary frequency spectra) 

Our aim is now to construct new analytical solutions to the Schroedinger equation, by 
following an exact (not approximate) approach. Let us, then, go back to eq.(l), and to 
its Bessel-beam solution (11), where, as before, relation (12b) holds: p p = ^2mE — p 2 z , 
with E = uh. 

The condition for obtaining a Localized Solution (cf. Fig. 3) is that 

E = Vp z + b , (33a) 

with b a positive constant (bearing the dimensions of an Energy, and regulating the 
position of the chosen straight-line in the plane (E,p z )); which corresponds in particular, 
on using eq.(12b), to the adoption of the integration limits 
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Figure 3: This time, the parabola and the chosen straight-line have equations E = 
p 2 z /(2m) and E = Vp z + b, respectively. The intersections of this straight-line with 
the parabola are now two, whose corresponding values are given in eq.(33b). Inside the 
parabola p 2 > 0. 

E ± = mV 2 (l ± yjl + ^s) + b . (33b) 

Localized Solutions can therefore be obtained by the following superpositions (inte- 
grations over the frequency, or the energy) of Bessel-beam solutions: 

z, = ew z f E+ dE Upvp/Ti) S(E) e*&< , (34) 

together with 

p p = ^ ^-E 2 + (2mV 2 + 2b)E - b 2 . (35) 

Notice that the in eq.(34) [as well as in eq.(39) below], the solution \l> depends on z, 
besides via (, only via a phase factor; the modulus of ^ goes on depending on z (and 
on t) only through the variable £ = z — Vt. 

4.1 Particular exact Localized Solutions 

We want now to re- write the integral X appearing in the r.h.s. of eq.(34) so that its 
integration limits are —1 and +1, respectively; that is, in the form 
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X = f 1 du S(u) Jo(^^VT^) e i/(c)u , 



quantity f(() being an arbitrary dimensionless function. To obtain this, we have to look 
for a transformation of variables [with A and B constants, with the dimensions of an 
Energy, to be determined] 

E = Au + B (36) 

such that 

p 2 p = P(l-u 2 ); u + = l; u_ = -l, (36') 

P being a suitable constant (with the dimensions of an Impulse square). On writing 
V 2 p 2 p = E(hV 2 M-E)-b 2 , with KM = 2m + 2b/V 2 , after some algebra one finds that 
it must be 

A = VPV; B = mV 2 + b; P = m 2 V 2 + 2mb. (37) 

Indeed, one can verify (by some more algebra) that eqs.(36)-(37) imply, as desired, that 

U- — — 1 and u + = 1. 

In conclusion, the transformation 



/ 2b 

E = mV 2 \ 1 + u + mV 2 + b (38) 

V mV 2 

does actually allow writing solution (34) in the form [recall that E = Au + B dE = 
Adu] 

tf(p,77,C) = M AeT' j l du S(u) Jo^vWT 3 ^) e^ u , (39) 



with 

T) = Z — Vt , 



16 



where v = V + b/(mV). Equation (39) is exactly, analytically integrable when S is a 
constant or a suitable exponential. 

Let us choose the complex exponential function (which will easily enter as an element 
in a Fourier expansion) 

S(E) = a n e^ nE , (40) 

with n an integer, and D = E + - E_ = 2mV 2 y/l + 2b/{mV 2 ), while Clfi £1X6 constant 
quantities (with dimensions of the inverse of an Energy). On remembering that E = 
Au + B, such a spectrum can be written in terms of u as 

S(u) = a n e iwnu e^ nB (40') 

(still with the dimensions of an inverse Energy). After some more algebra, the analytic 
exact solution to the Schroedinger equation, corresponding to spectrum (41), results to 
be[ll] 

*(p, n, C) = ATa n 2A — e l ° nB , (41) 

where A, B, P are given by eqs.(37) and 

Equation (41), as we have just seen, is a particular exact Localized Solution to the 
Schroedinger equation; but we are going to utilize it essentially as an element of suitable 
superpositions. Before going on, however, we wish to depict in Figs|4] an elementary 
solution: namely, the square magnitude of the simple solution corresponding, in eq.(34), 
to the real exponential 

S(E) = s exp[a(E-E + )] , (43) 

a being a positive number, endowed with the dimensions of an inverse Energy, as well as 
So- When a = 0, one ends up with a solutions similar to Mckinnon's[23]. Spectrum (43) 
is exponentially concentrated in the proximity of E + , where it reaches its maximum value; 
and becomes more and more concentrated (on the left of E + , of course) as the arbitrarily 
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chosen value of a increases. To perform the integration in eq.(34), it is once more useful 
to operate the variable transformation (36) and go on to eq.(39), spectrum (43) assuming 
now the form 

S(u) = s e~ aE+ e aB e aAu . 

Performing the integration in eq.(39), by a process similar to the one which led us to 
eq.(41), in the present case we get 

#(p,77,C) = VVs 2V/Pexp[z"^r/] exp[-aVVP] S ^ (44a) 

where 

Y = & Jpi-thaV + iO* , (44b) 

quantity P having been defined in eq.(37); and one should remember that rj = z — vt is 
a function of b. 

Equations (44) appear to be the simplest closed-form solutions (see Figsji]) to the 
Schroedinger equation, since they do not need any recourse to series expansions of the 
type exploited in the following Subsection. However, the solutions that we shell construct 
below can correspond to spectra more general than (43); for instance, to the gaussian 
spectrum, which possesses two advantage w.r.t. spectrum (43): it can be easily centered 
around any value of u, that is, around any value E of E in the interval [E_,E + ], and, 
when increasing its concentration in the surrounding of E, its "spot" transverse width 
does not increase, at variance with what happens for spectrum (43). Anyway, the exact 
solutions (44) are noticeable, since they are really the simplest ones. 

Some physical (interesting) comments on the results in eqs.(44) and Figs. 4 will appear 
elsewhere. Here, let us add only a few further Figures and some brief comments. Let us 
first recall that, as predicted in the first one of Refs.[6], the Localized (Nondiffracting) 
Solutions to the ordinary wave equations resulted to be roughly ball-like when their peak- 
velocity is subluminalfll], and X-shaped[4,7] when superluminal. 
Now, normalizing p and (, we can write eq.(44b) as 
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(a) (b) 

Figure 4: In these figures we depict an elementary solution: namely, the square magnitude 
of the simplified solution, eq.(44a), corresponding to the real spectrum S(u) = s exp[(E — 
E + )a\i as a function of p' = py/ P '/% and of (' = (y/P/h. Quantity a is a positive number 
[when a = one ends up with a solutions similar to Mckinnon's[23]], while b for simplicity 
has be chosen equal to zero. Figure (a) corresponds to a = E+/5, while figure (b) 
corresponds to a = 5E + . For the properties of the spectral function (43), see the text. 

with p' = y/Pp/h and (' = \/P(/h, quantity P being given by the last one of eqs. (37), 
namely P = m 2 V 2 + 2mb, while A = aA = y/PaV. For simplicity, let us confine 
ourselves to the case 6 = 0, forgetting now about the more interesting cases with 6^0; 
therefore, it will hold the simple relation 

A = maV 2 . 

In the present case of the Schroedinger equation, we can observe the following. 

If we choose A = 0, which can be associated with V = 0, we get the solutions in 
Figs. 5: that is, a ball-like structure. 

By contrast, if we increase the value of A, by choosing e.g. A = 20 (which can 
be associated with larger speeds), one notices that also a X-shaped structure starts to 
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M 2 




Figure 5: In these, and the following Figures 6, 7 and 8, we depict the square magnitude of 
some more solutions of the type (44a), normalized with respect to p and £; still assuming 
for simplicity b = 0, so that A = maV 2 . The present figures show the "ball-like" structure 
that one gets, as expected, when A = (see the text, also for the definitions of p' and ('). 
Fig.(b) shows the projection on the plane (//, (') of the 3D plot shown in Fig. (a). 

contribute: See, e.g., Fig. 6. 




Figure 6: The solution, under all the previous conditions, with an increased value of A, 
namely with A = 20. An X-shaped structure starts to appear, contributing to the general 
form of the solution (see the text). 
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To have a preliminary idea of the "internal structure" of our soliton-like solutions to 
the (ordinary) Schroedinger equation, let us plot, instead of the square magnitude of 
its real or imaginary part: Let us choose its real part, or rather the square of its real 
part. Then even in the A = case one starts to see the appearance of the X shape, 
which becomes more and more evident as the value of A increases: In Figs. 7 we show the 
projections on the plane ((', p') of the real-part square for the solutions with A = 5 and 
A = 50, respectively. Further attention to such aspects will be paid elsewhere. 



Re(T) Re(^) 




(a) (b) 

Figure 7: To get a preliminary idea of the "internal structure" of our soliton-like solutions, 
it is useful to have recourse (see the text) to the real part of ^. In these Figures we plot 
the projections on the plane ((', p') of the real-part square for the solutions with A = 5 
(figure (a)) and A = 50 (figure (b)), respectively. 



But the (square of the) real part of \1/ does show, in 3D, also some "internal oscil- 
lations": Cf., e.g., Fig.8 corresponding to the value A = 5. We shall face elsewhere, 
however, topics like their possible connections with the de Broglie picture of quantum 
particles, et alia. 
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Figure 8: The (square of the) real part of \I> shows, in 3D, also some "internal oscillations": 
this Figure corresponds, e.g., to the value A = 5. 



4.2 A general exact Localized Solution 

Let us go back to our spectrum S(E) in eq.(40). Since in our fundamental equation 
(34) the integration interval is limited [E- < E < E + ], in such an interval any spectral 
function S(E) whatever can be expanded into the Fourier series 

00 

S(E) = £ a n e^ nE , (45) 

n=— oo 

with 

G - = ^f E+ dES(E)e~^ nE , (46) 

JJ JE- 

quantity S(E) being an arbitrary function, and D being still defined as D = E + — E-. 

Inserting eq.(45) into eq.(34), and following the same procedure exploited in the 
previous Subsection (in particular, going on again from E to the new variable u), we end 
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up — after normalization — with the general exact localized solution to the Schroedinger 
equation: 

*(p,r],0 = M2Ae^ £ a n exp [i^nB] ^ , (47) 

n=-oo U Z 

where Z is defined in eq.(42), and the coefficients a n are given by eq.(46). 

It is worthwhile to note that, even when truncating the series in eq.(47) at a certain 
value n = N, the solutions obtained is still an exact LS of the Schroedinger equation! 



5 About finite-energy Localized Solutions to the 
Schroedinger equation 

The solutions found above, even if very instructive, are ideal solutions which are not 
square integrable; and cannot be accepted in QM. It is important, therefore, to show how 
to construct finite-energy solutions. 

Let us obtain localized solution to the Schroedinger equation endowed with finite 
energy, by starting from eqs.(44). First of all, one has to integrate over b by adopting 
a spectrum S(b) strongly bumped around a value b : We already know, indeed, that 
spectra of this type are required in order to get solutions that are non-diffracting all along 
a certain field-depth. 

Then, it can be easily seen that the finite-energy solution, \l/f c , can be preliminarily 
written as 

*fe = M (/_ - I + ) , (48) 

where J_ and I + are two (dimensionless) integrations over b from to infinity (quantity b 
having been defined in eq.(33a), and therefore having the dimensions of an Energy), while 
s appears in eq.(43). 

Let us now pass from b, defined in eq.(33a), to the new variable w = \[P . One has to 
choose a spectrum S(w) corresponding to a S(b) concentrated around a specific value of 
b; let us therefore adopt the gaussian function 
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S(w) 



m^fq 



exp[— q(w — w o y 



(49) 



with w > mV > 0. 

When we go on from b to the new variable w = y/P (where P depends on b), the 
two quantities i_ and I + become integrations over w from mV to oo. After further 
calculations, and using relation 3.322.1 in ref.[22], one obtains that 



I± = 



u 



-qwo 



e 2R exp [ 



wl 

JJ2 J 



(50) 



where 



U = 2\q + 1; W± = -2qw + aV±i^ , 

V 2m y/P 

quantity Y having been defined in eq.(44b). 

We have therefore shown that realistic (finite-energy) Localized Solutions exist also 
to the Schroedinger equation; they will be non-diffracting only till a certain finite distance 
(depth of field). The analysis of explicit, particular examples will be presented elsewhere. 



6 The case of non-free particles 

Let us consider now the case of a particle in the presence of a potential: for simplicity, 
let us confine ourselves to the case of a cylindrical potential. 

Namely, let us consider the Schroedinger equation with a potential of the type U(p): 

Now, we can use the method of separation of variables writing ip = R(x,y)Z(z)T(t). 
With this, we get the well known solutions 
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T = e-* Et 



(52) 



2 = e Wzz/h 



(53) 



and the eigenvalue equation 



h 2 W\R + 2mU(p)R = A 2 R 



with 



A 2 = 2mE-p 2 z 



(54) 



(55) 



Supposing a potential U(p) that only allows transverse bound states (as the parabolic 
potential), we will find eigenfunctions R n (x,y) and discrete (degenerate) eigenvalues A 2 . 
We can construct more general solutions 



m = Y,fnRn{^yy kzZ,n e^ Et 



(56) 



with 



2mE = p 2 z + A 2 n 



(57) 



Considering p z > (forward propagation), the constraint (57) defines a set of parabo- 
las (something like the modes in a waveguide: Cf. Refs.16). Chosen a certain A 2 , once a 
value for p z is given, the value of E gets fixed. 



To obtain from (56) a train of localized pulses, i.e., a wavefunction ty(x,y,z — Vt), 
we must have 



E = Vp z 



(58) 



So, from conditions (57) and (58), p z must assume the values 



p z = mV \l±\ 1- 



m 2 V 



A 2 , 



2T/2 « 



(59) 
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with 



A n < mV (60) 
Figure 9 illustrates the situation. The values to E and p z that furnish localized pulse 



trains are given by the intersection between the parabolas defined by eq.(57) and the 



straight line defined by eqs.(p8|. Note that in these cases the series (56) will be always 



truncated (finite number of terms), due the condition (60). We also have to note that, for 



any given A^, one gets two possible values of k z (see eq.(59)), as it can be observed from 
Fig.9, in which the straight line cuts each parabola twice. 




Figure 9: In the case of a particle in the presence of a cylindrical potential, the values to 
E and p z that furnish Localized Pulse trains are given by the intersection between the 



parabolas in eq.(57) and the straight line in eq.(58): see the text. It can be noticed that, 



for any given A;;, one gets two possible values of k z (cf. eq.(59)), since the straight line 
cuts each parabola twice. See the text, and cf. also Refs.[17]. 



For our purpose, the superposition has to be 



V(x,y,z-Vt) = J2fn R n(x,y)e 



ip z „(z-Vt)/h 



(61) 



with 
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Pz = mV|l±Wl-^A*| (62) 



and 



A n < mV (63) 

In principle, any set of coefficients /„ will furnish trains of localized waves. 
Observation!: If we look for a square- integrable wave function, we can start from 



superposition (56) and integrate its terms over p z around each p zn , respectively (as we 
already did in our papers on X-type pulses propagating along wave-guides[17]). But in 
the present case, in general, the group-velocities defined at the points p zn will not be 
the same, as it happened in the waveguide case; and we will therefore meet a kind of 
intermodal dispersion, besides the group- velocity dispersion. Let us recall, incidentally, 
that such an intermodal dispersion did not occur in the case of X-type waves, traveling 
in metallic wave-guides, due the peculiar fact that the group-velocities defined at those 
points were always the same ). After the integration, we can obtain an envelope with a 
train of pulses (or just one pulse) inside it. The envelope will suffer dispersion, but the 
train of pulses inside it will not. 

More general localized wave trains can be obtained using the relation E = Vp z + b, 
with b a positive constant. 

In the case of potentials like U(p), one can search for solutions with cylindrical sym- 
metry, for simplicity. However, solutions without this symmetry can be investigated: and 
they will be interesting for an analysis of angular momentum. 
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